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Abstract - Current models of the effect of spontaneous emission on the electron beam dynamics 
neglect the discreteness of electron recoil associated with photon emission. We present a novel, 
one-dimensional model of the effect of spontaneous emission on the electron beam dynamics in an 
undulator both in the classical regime where discrete electron recoil is negligible, and the quan- 
tum regime where it is significant. It is shown that in the classical regime, continuous decrease 
of the average electron energy and diffusive growth of the electron energy spread occurs, in agree- 
ment with previous classical models. In the quantum regime, it is shown that the evolution of 
the electron momentum distribution occurs as discrete momentum groups according to a Poisson 
distribution. The narrow momentum features of the quantum regime may be useful for generation 
of coherent radiation, which relies on electron beams having sufficiently narrow momentum/energy 
distributions. 



1. Introduction. — Free electron lasers (FELs) using relativistic electron beams in 
static undulators or wigglers have recently produced high-intensity coherent radiation in 
the XUV and X-ray regions of the spectrum [TJ[2] . Methods for extending the wavelength 
range of coherent FEL-like sources to even shorter wavelengths using magnctostatic or elec- 
tromagnetic laser wigglers arc being actively investigated [3J- These efforts require us to 
revisit theoretical models used to describe the emission of radiation by relativistic electrons. 
In most current models of coherent radiation generation by e.g. the Free Electron Laser 
(FEL), the discreteness of the recoil due to spontaneously emitted photons can be safely 
neglected. In previous classical theories [U[5] the momentum spread produced by sponta- 
neous emission is described as a continuous diffusive process. At very short wavelengths 
(high photon momentum/energy) the discreteness of the momentum recoil will eventually 
become significant. The example we consider is that of a relativistic electron beam passing 
through a magnetic undulator but the analysis is also applicable to electromagnetic laser 
wigglers and betatron oscillations in plasma channels [5] as in all cases the electrons oscillate 
transversely and consequently radiate spontaneously. This paper describes a model which 
allows a description of the effect of spontaneous emission on the electron beam dynamics in 
both the classical regime where discrete electron recoil is negligible, and the quantum regime 
where it is significant. In section 2 we describe this model in the limit of weak wiggler fields 
(wigglcr parameter a M C 1) where emission occurs only at the fundamental. In section 3 we 
demonstrate the main features of the classical and quantum regimes of emission, and in sec- 
tion 4 we describe how the model can be extended to include emission at higher harmonics, 
which are significant for a w > 1. 
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2. Spontaneous Emission Model. — We use a model based on the evolution of a 
distribution function for the electron momenta. In this section, we consider spontaneous 
emission at the fundamental frequency alone, and harmonics are neglected. The effect of 
higher harmonics is described in section 4. As was shown in [7], in the limit where a w < 1, 

2rraa 2 

an oscillating electron spontaneously emits photons at a rate R = 3A m per unit distance 
through the undulator, where a is the fine structure constant, the wiggler 

parameter and A™ is the undulator/wiggler period. As the electron beam passes through 
the undulator, photons are emitted at random times. Consequently after a distance z, 
different electrons will have emitted different numbers of photons. If the beam is initially 
monoencrgctic there will therefore be a growth of energy/momentum spread i.e. so-called 
'quantum excitation' of the beam in addition to a decrease in the average beam momentum 
i.e. 'radiation damping'. The novel feature of our model is that it takes account of the fact 
that spontaneous emission involves emission of photons with momentum hk, where k = 2tt/X 
is the photon wavenumber and A is the wavelength of the emitted radiation which in the 
case of the FEL can be expressed as A = A.^ ^ ^-y™ ) • Consequently the probability of an 
electron having momentum p will be increased by spontaneous emission from electrons with 
momentum (p+hk) but decreased by spontaneous emission from electrons with momentum 
p. From this physical picture we can write a rate equation for the distribution function due 
to spontaneous emission alone: 

dw(z,p) =RW { ZiP + nk)-RW(z,p) (i) 

oz 

where W(p, z)dp is the probability of finding an electron between momenta p and p + dp. 
In order to see the statistical behaviour which this model describes, it is possible to rewrite 
the electron momentum in terms of number of photons emitted, N, where N is a statistical 
variable. In order to do this, we use the dimensionless momentum variable p' = p/hk = 
H — N. Consequently, eq. ([I]) can be rewritten as 

= RW(z,N -l)-RW(z,N) (2) 

which is the master equation for a Poisson process |S]. This has a solution 

< N > N e- <N> 
W(z,N) = ^ W] 

where the mean and variance of N are described by 

(N) = ((N - (N)f) = Rz 
. Consequently the electron momentum p = po — Nhk will behave as 

(p) = po — hkRz 

((j>-(p)f) = {hkfRz 

where for simplicity we have assumed a beam which is initially monoenergetic. It should be 
noted that our model neglects the finite linewidth of the spontaneously emitted radiation 
i.e. it assumes emission of photons of a well-defined momentum hk. As the linewidth of 
wiggler radiation is ^j- ~ , then the corresponding uncertainty in photon momentum is 

hAk ~ jj*- . Consequently the neglect of the finite linewidth of the radiation in our model 
is valid if N w 3> 1. For consistency, this means that the probability of photon emission per 
wiggler period is small i.e. RX W ~ 2aa^ « 1, which is satisfied if <§C ~ 70. This 
assumption is therefore valid for most FELs which use wigglers with a w ~ 1. 
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3. Classical and Quantum Regimes. — In order to observe the different regimes of 
spontaneous emission, we use the dimensionless variables P = — , e = — and Z = Rz, where 

(7 p (Tp 

a p is a characteristic momentum spread of the system e.g. that produced by interaction or 
instability such as in the Free Electron Laser (FEL). Note that the characteristic momentum 
spread, <r p , in the case of a low-gain FEL would be that associated with the lincwidth of 
the spontaneous emission curve i.e. a v = a 1 mc = -£p£, where m is the electron rest mass, 
c is the speed of light and N w is the number of undulator/wiggler periods. In the case 
of a high-gain FEL, a p = a y mc = pr/mc, where p is the classical FEL parameter and the 
parameter e becomes 

1 

P 

where p = ^j^- p is the quantum FEL parameter [9l[T0] . 
In terms of these variables eq. [T] becomes 

^fz P) =W(Z,P + e)-W(Z,P). (3) 

According to the analysis of cq. [5] above, the average value of P will decrease linearly at rate 
e i.e. 

(P) = (P ) - eZ (4) 

and P will diffuse according to 

a 2 = {{P-(P)f) = e 2 Z (5) 



3.1 The Classical Limit (e — > 0). In the limit e — > it is possible to approximate the 
finite difference equation in eq. [5] using a Taylor expansion and retain only the first two 
terms i.e. 

^,P +e ,^ (z ,P )+e ^4«2 + ... 

Consequently, the finite difference equation in cq. (|2|) reduces to the drift-diffusion equation 

dWjZ, P) _ dW(Z, P) e 2 d 2 W(Z, P) 

dZ ~ £ dP + ~2 W 2 ■ () 

Converting to unsealed variables, the expression in eq. [S]can be written as 

^ A72) ^0 2i?z = 3i fc ^ a » 74z (7) 

which is approximately the same expression (to within a numerical factor ~ 1) as those 
derived in [4] and [5] when a w < 1. 

A numerical example of classical evolution is shown in figure [1] for the case where e = 0.2. 
The drift and diffusive growth towards a Gaussian-like momentum distribution can be seen 
clearly as Z increases. 

3.2 The Quantum Limit (e > \). In the case where e is sufficiently large that the 
higher-order terms in the Taylor expansion of W(Z, P + e) must be retained, the reduction 
to a drift-diffusion equation can no longer be performed and the finite difference form of 
cq. ([2]) must be retained. A numerical example of quantum evolution is shown in figure [2] 
for the case where e = 10 using the same initial condition as that used for fig. [TJ It can be 
seen that in contrast to the classical case discussed previously, the momentum distribution 
consists of a series of discrete lines separated by e (hk in unsealed units). It should be noted 
that while the envelope of the distribution broadens with increasing Z according to eq. (0]) 
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and |[SJ| , the width of the individual discrete momentum groups do not, remaining close to 
their original values. 

Note that even for e > 1, in order to observe the quantum discreteness of the momentum 
distribution, it is necessary that the initial momentum spread , <7o, be sufficiently small that 
the spacing between the momentum lines can be resolved i.e. ctq < e. This is confirmed in 
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Fig. 3: The electron momentum distribution at Z — 5 and the evolution of the mean and variance 
of the momentum distribution as calculated from a numerical solution of eq. [5] when e = 0.2 and 
(70 = 10. 



4. Spontaneous Emission Including Harmonics. — In the case of a planar undu- 
lator where the undulator parameter a w > 1, the emission of radiation at odd harmonics of 
the fundamental frequency can no longer be neglected. We now extend the model derived 
in section n to include the effect of radiation at higher harmonics. 

When electrons radiate higher harmonics in addition to the fundamental, then the 
spontaneous emission process involves emission of photons with momentum nhk, where 
n = 1, 3, 5, ... Consequently the probability of an electron having momentum p will be in- 
creased by spontaneous emission from electrons with momentum (p + nhk) but decreased 
by spontaneous emission by electrons with momentum p. From this physical picture we can 
write a rate equation for the distribution function due to spontaneous emission including 
higher harmonics analogous to eq. ([1} : 



dW(z,p) 



n=l,3,5, 



[R n W(z,p + nhk) - R n W(z,p)} 



[RnW(z,p + nhk)}-W(z, P )i J2 Rn ) ( 8 ) 

n=l,3,5,... \n=l,3,5,... / 

where R n is the spontaneous emission rate per unit distance of electrons with momentum 
nhk. R± is the fundamental spontaneous emission rate and is related to the rate defined in 
section n by R\ — j^t (Jo(C) ~ Ji(C)) 2 where ( = 2 {i+a' i ) so R\ ^ R when o„«l 
[llj . The emission rate at harmonic n is 

p»dM- J?±±(n()\ 2 

Rn -{ 2 MO -MO ) Rl - 
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Repeating the procedure of section „ it is possible to rewrite eq. ([8]) in terms of the dimen- 
sionless variables Z = R\Z, P = and e = where it should be noted that Z and P 
have been scaled with respect to the fundamental emission rate (i?i) and photon momentum 
(hk) respectively. Consequently the rate equation including harmonics can be written as 



dW(Z, P) 
dZ 



E 

n=l,3,5,. 



[Cl n W(Z ) P + ne)-fl n W(Z,P)} 



p n W{Z,P + ne)]-W{Z,P){ ^ 



n, 



(9) 



n=l,3,5,... 



^n=l,3,5, 



where 



n = — 



Pi, 
Pi 



'Jn-iinC) - Jn±i(n<;y 

2 mo- MO 



Fig. |4] and fig. [5] show the influence of harmonics on the evolution of the electron momentum 
distribution for different wiggler parameters, a w , when e = 0.2 and e = 10 respectively. 
Fig. [5] and fig. [7] show corresponding evolution of the mean momentum and the momentum 
variance for different wiggler parameters, a w , when e = 0.2 and e = 10 respectively. It 
can be seen from these figures that in both classical and quantum regimes as a w increases, 
both the rate of decrease of the mean momentum (drift) and the momentum variance also 
increase. The reason for this is that when the electrons emit a harmonic, n (where n > 1), 
each spontaneous emission event causes population exchange between momentum groups 
separated by nhk, which are more widely separated in momentum than in the case of 
emission at the fundamental alone (where n=l). 
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Fig. 4: Classical regime including emission at harmonics : Momentum distribution at Z — 10 
calculated from a numerical solution of eq. @ including harmonics when e = 0.2 and oo =0.1 for 
various wiggler parameters a w : (a) a w = 0.1, (b) a w = 1 & (c) a w = 2. 



5. Conclusions. — We have presented a novel dynamical model of spontaneous emis- 
sion by a relativistic electron beam spontaneously emitting undulator radiation, which is 
applicable to both magnetostatic wigglers and electromagetic laser wigglers (i.e. Thomson or 
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Fig. 5: Quantum regime including emission at harmonics : Momentum distribution at Z — 5 
calculated from a numerical solution of eq. including harmonics when e = 10 and <jq = 0.1 for 
various wiggler parameters a w : (a) a w — 0.1, (b) a w = 1 & (c) a w — 2. 



Compton scattering) in both classical and quantum regimes of emission. Using this model, 
we have shown that in the classical limit where discrete photon recoil is negligible, the model 
predicts a linear drift and diffusion of the momentum distribution at a rate which agrees 
well with those of previous work using alternative classical models [H[5]. I n addition, our 




Fig. 6: Evolution of mean momentum and momentum variance when e — 0.2 and <to = 0.1 (i.e. 
classical regime) for various wiggler parameters. 
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Fig. 7: Evolution of mean momentum and momentum variance when e = 10 and oo =0.1 (i.e. 
quantum regime) for various wiggler parameters. 



model shows that in the quantum limit where photon recoil is significant the evolution of the 
momentum distribution changes dramatically to involve discrete narrow momentum groups 
under an envelope which drifts and diffuses. When emission at harmonics is also significant, 
it was shown that as a w increases, both the rate of decrease of the mean momentum (drift) 
and the momentum variance also increase. 

To our knowledge this is the first time that the dynamical evolution of electron beam 
dynamics due to quantum effects arising from spontaneous emission by a beam of oscillating 
electrons have been consistently described. By combining this model with the coherent emis- 
sion processes in e.g. an FEL we should now be able to get a truer picture of the significance 
of quantum effects associated with spontaneous emission. Of particular interest is its impact 
on coherent FEL emission at short wavelengths using classical and quantum modes of FEL 
operation. To date these considerations have been based on the evolution of the envelope of 
the momentum distribution, which in the quantum regime can completely neglect the much 
narrower momentum features associated with discrete photon recoil. As coherent radiation 
generation typically relies on electron beams with narrow momentum spreads, this difference 
could be significant. The model described here retains the narrow features of the quantum 
regime, so consequently should be useful for investigation of schemes in which the aim is 
control and exploitation of the discrete structure of the electron distribution arising from 
the quantum regime of spontaneous emission for coherent short-wavelength generation. 
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